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1 Introduction 



There is now strong evidence that string theory contains a large number of metastable vacua, 
collectively referred to as the landscape [1]. In order to make progress in understanding the 
place of a typical observer in this landscape of possibilities we must understand which vacua 
are possible, and the method of their population. The former of these problems is a current 
topic of research and we will not attempt to contribute to this field here. Instead we will take 
the pioneering model of Bousso and Polchinski [2] , and investigate the latter. 

The population mechanism for the landscape is eternal inflation. De Sitter vacua are con- 
stantly expanding and, in string theory, are meta-stable. As these vacua expand they nucleate 
bubbles of new vacuum, which themselves expand and nucleate further bubbles. In this way 
it is expected that all corners of the landscape are realiseci^, and the pertinent question is 
then which are the most probable. The fundamentals involved in answering this question are 
currently a matter of some contention: A difficulty often referred to as 'the measure problem.' 

Various different measures have been proposed in the literature. We will not attempt to 
describe them all here, nor compare their various implicit benefits and difficulties. Instead, 
we will take one particular class of measures, and investigate the outcome of using them to 
predict the likely vacua that observers will find themselves in. The class we will study are the 
volume weighted measures, proposed by Linde [1]. This class embodies the appealing intuition 
of rewarding the fastest expanding vacua. A refinement of this measure has recently been shown 
to produce results with surprisingly little dependence on the choice of time parameterization, 
and that do not exhibit a hotness (or youngness) paradox [5j. Using the Bousso- Polchinski 
landscape to investigate the results of volume-weighted eternal infiation has the advantage of 
this model already having been used to investigate other measures, namely the bubble-counting 
measure of [HI E] and the holographic measure of [HI E] • This will allow us to directly compare 
the results of several candidate measures within the same theoretical framework. 

The Bousso-Polchinski landscape has inherent in it a number of approximations. For ex- 
ample, in this toy model the charges and their associated moduli are independent of fiuxes, 
which is unlikely to be true of the real landscape. The contribution of slow-roll infiation is also 
ignored. Within the Bousso-Polchinski framework we will make several further approximations, 
for analytic convenience, such as treating discrete fiuxes as continuous variables. Nonetheless 
we will be able to find, within this model, reliable predictions for the distribution of the cos- 
mological constant, the likely configuration of fiuxes and even the typical lifetime of a 'pocket 
universe'. 

In section 2 we will give a brief account of volume weighted eternal infiation and, using 
examples, will give an account of how to calculate probabilities in this measure. We will follow 
this with a description of the Bousso-Polchinski landscape, and the transition rates expected 
within it. We will then refiect on some anthropic considerations that can be applied to this 
landscape. In particular, we consider anthropic bounds on the cosmological constant and on 

'''For a recent study of possible mechanisms to divide the landscape into different regions see [3]. 
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the lifetime of vacua that must be obeyed if observers are to be expected to occur. 

Section 3 gives a description of how the high energy vacua in this landscape are expected to 
be populated in volume weighted inflation. We discuss here how initial conditions are quickly 
forgotten, once the highest energy, most rapidly expanding vacua are realised, and subsequently 
come to dominate the evolution of the universe. We then argue that above a certain high 
energy boundary the distribution of probability is approximately constant, and isotropically 
distributed. 

In section 4 we consider tunneling events down from this high-energy barrier to a surface one 
jump up from the anthropic shell. We find that the primary contribution to probability on this 
low energy shell is due to chains of transitions occurring in a straight line, in one dimension 
of the fiux space. Approximating the discrete space as continuous, we then find an explicit 
expression for the total tunneling rate down from high energies to the penultimate sphere, 
before the anthropic one. We find that there exist peaks of probability on this penultimate 
sphere, and find their position at various different levels of approximation. 

Section 5 gives a description of the final transition down into the anthropic shell. This 
transition is treated separately as the effects of gravity are negligible for this flnal jump. We 
find that this last step also produces peaking, but that this peaking is negligible in comparison 
to that found in section 4. 

A description of the extent of peaking, due to suppressed transitions, can be found in section 
6. For a typical set of parameters we find that 99.9999% of the vacua in the anthropic shell 
should be expected to lie within a fraction ~ 10~^^ of its surface area. The selected vacua have 
one fiux much larger than the others. There are ample vacua in the selected set to allow for an 
anthropic solution of the cosmological constant problem. We also consider in this section the 
effect of the anthropic life-time bound, already discussed in section 2. The peak of probability 
is found to be centered on vacua with insufficiently long life-times for observers to appear. If 
we consider only vacua with life-times > 10 billion years, then the peaking in probability exerts 
a pressure favoring the vacua with shorter lives. We find that a typical vacuum, in which 
observers take 10 billion years to appear, will be most likely to decay approximately 2.2 billion 
years later. 

Section 7 gives a discussion of the results obtained, and how they may be altered by effects 
that were not included in our model. For example, if multiple branes are allowed to be nucleated 
simultaneously then the peaking we have found should be expected to be reduced, and to vanish 
in the limit of indefinitely many branes nucleating simultaneously. We also consider how the 
required stage of slow-roll infiation, before re-heating, could modify our results. 



2 The Model 

We consider here volume-weighted measures of inflation in the Bousso-Polchinski landscape. 
We will begin with a brief explanation of these measures in eternal inflation, and then describe 
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Figure 1: A toy landscape, with three dS vacua. Kij denotes the tunnehng rate from vacuum i 
to vacuum j. 

the features of the Bousso-Polchinski landscape that are of relevance for our study. 
2.1 Volume Weighting 

The essential feature of volume weighted measures is that they reward those regions which 
expand most rapidly. That is, in a given increment of time. At, the 3-volume of each type of 
vacuum is increased as a function of its rate of expansion, in some proportionate way. Clearly 
the way in which time is measured is of primary importance here. For definiteness we will define 
'time' as being the proper time measured by a congruence of comoving observers. This can 
be measured in an invariant way in units of the 11-dimensional Planck time, or by equipping 
each observer with a suitably constructed light clock. Inflation occurs in each vacuum until a 
tunneling event creates a bubble of new vacuum in the old (or on rare occasions until Minkowski 
space is reached, if this exists in the theory). The number of points in any particular vacuum 
that do not eventually undergo a transition to a bubble of new vacuum are a set of measure 
zero. However, due to the rapid expansion of high-energy vacua, such vanishingly small sets 
can still come to dominate the evolution of the universe. To make these ideas concrete, we will 
illustrate this behavior with an example. 

Consider a simple landscape of 3 de Sitter vacua, as shown in Figure [TJ Here denotes 
the transition rate from vacuum i to vacuum j, and the Hi are Hubble constants. We assume 
that vacuum 1 has a much higher cosmological constant than vacua 2 and 3, so that Hi ^ if2,3 
and Kii ^ Kii, K,23, 1^32- If Pi IS the (unnormalized) probability of finding a given point in the 
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vacuum state i, then we can describe this dynamical system through a set of coupled ODEs 



Pl = f^2lP2 + K31P3 - (/tl2 + Ki3)Pl + SHiPi , 
P2 = K12P1 + KS2P3 - (/t23 + /t2l)^2 + 3i72P2 , 
h = f^ViPl + /«23^2 - (/«32 + taO^S + 3i^3P3. 

As long as Hi ^ for all i, j (i.e. the highest energy vacuum expands much faster than it 
decays), the first of these equations becomes 

Pl - 3HiPi Pl oc e-^^'K 

The terms dominating the RHS of the second two equations are then those which are propor- 
tional to Pl, as t ^ 00. These equations are solved by 

Pi f^li -. Pi f^li 

— ~ ^1 or ~ . 

Pl 3Hi Pj Kij 

This confirms that the probability of finding a point in any particular vacuum is proportional 
to its rate of decay from the highest energy (fastest-growing) vacuum. Essentially, the most 
rapidly expanding vacuum sources all of the others - it is easier to make the lower vacua by 
expanding in vacuum 1 and then decaying down, than it is to do so by simply expanding with 
the smaller Hubble constant of the lower vacua. 

These results readily generalize to more complicated systems. First, consider a landscape 
in which the highest energy vacuum can decay to all other vacua. In this case we have 

Pj ^hi -, Pi ^hi /-, \ 

7r-^<l or 7^-— ' (1) 

-T h -J-n/i JTj Kfij 

where h labels the highest energy vacuum. Now consider a landscape in which there exist 
vacua that cannot be created from direct decays from the highest energy state. In this case the 
probability of finding a point in the vacuum state j is 

P, ^ 3Hk 3H, 3H, 3H, ^ ^ 

all 

paths 



h „ J-n/i onh orih 

all 



The chain of K,ij/3Hh gives the probability of starting in the state h and finishing in j, following 
some series of decays through the vacua a, b, c... To obtain the total probability of finding a 
point in the vacuum j it is then necessary to sum over all possible paths. Clearly some paths 
will provide more of a contribution than others to this sum. In fact, as long as Hij <^ Hh, 
only the shortest decay chains will give a significant contribution. Every chain with additional, 
unnecessary, intermediary vacua will be suppressed by factors of Kij/3Hh -C 1. However, if 
Kij ~ 3Hh then the products of decay are no longer suppressed. In this case the progeny of 
decaying vacua will be as equally abundant as the vacua from which they decayed. Rapid 
transitions of this type will not contribute suppressing factors to the decay chains mentioned 
above. 
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2.2 The Bousso-Polchinski Landscape and Transitions 



We will consider here the Bousso-Polchinski (BP) landscape of M-Theory [2] . This model has 
the advantages of being relatively simple and, most importantly for our purposes, having known 
expressions for the transition rates. 

The BP landscape consists of a discrete set of vacua, each with a different effective cosmo- 
logical constant given by 

1 ^ 

^ = -^o + ^Y.^lnl (3) 

i=l 

where J is the number of fluxes, qi the charge associated with each of them and the number 
of units of a particular flux. The discreteness of the rii arises from the quantization of four-form 
fluxes in M-theory compactifications, and their multiplicity from the geometry of the M-theory 
compactification - each non-trivial three-cycle giving rise to a different four-form flux. The 
number of such cycles is generically quite large. This multiplicity allows a much finer spacing 
of A than that given by a single flux alone, which typically has a charge only a few orders of 
magnitude below the Planck scale [2]. 

From equation ([3]) it is clear that for a given choice of A, the fluxes are constrained to lie on 
a sphere in flux space, of radius a/2(A + Aq). Taking ^ A < 10~^^° (in 4D Planck units) we 
obtain a thin shell of radius ~ \/2Ao, which corresponds to the set of vacua in the BP landscape 
which have an anthropically acceptable cosmological constant (see below). 

Transitions in this landscape occur via the nucleation of branes. The transition rates, 
including gravitational effects, are given by [12] 

= AijC-^^^ . (4) 

The Aij can be approximated by Aij ~ 1 (any effect they have will be exponentially unim- 
portant). Now let us consider the nucleation of bubbles which change the flux through some 
particular cycle by one unit, from rij to nj — 1, such that 

|AA,| = {n, - l/2)g2 . 

The tunneling rate is then given as the product of its value in fiat space, times some gravitational 
correction, as [13] 

i?,,=i?J*^P-V(x,y), (5) 
where the gravitational correction factor is 



and the fiat space rate is 



pflat-space 27?! 1 

""^i 8 (n.-l/2)3gf 
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The dimensionless parameters in r{x,y) are given by 



X 



8IAA,- 



K^* - 1/2) 



and 



y 



2A 
IAA" 



For our purposes rii will always be ^ 1, so that x ~ 3/8n -C 1. The variable y has a clear 
interpretation as the energy-level of the parent vacuum. At high energies Aj ^ |AAj|, so that 
y ^ 1. Conversely, at the lowest energies Aj ~ |AAj|, so that y ~ 1. The combination xy can 
be interpreted as a measure of distance from the axes of the flux space. Consider the surface 
defined by x?/ = 3/8, we then have 



xy 



3 2A, 
8ra,- IAA,' 



3 A 



Anfq'^ 8 



E2 2 



2An 



The LHS of this second equation is the perpendicular distance (squared) away from the axis 
corresponding to the direction in which the transition is occurring. The RHS is the radius 
(squared) of the anthropic shell. Therefore when xy < 3/8 the transition is occurring in 
a cylinder around this axis, of the same cross-sectional area as the anthropic shell. When 
xy ^ 3/8 the transition is occurring well away from the axis. 

Useful estimates of tunneling rates can be found in a number of different regimes [7]. Firstly, 
at high energies far from the axis, when y ^ 1 and ^ 1, it can be seen that 



r{x,y) ~ V2{xy) < 1. 



(9) 



The tunneling rate ([5]) is then given by 



Kij ~ Exp 



A 



3/2 



(10) 



This expression is independent of rij. Alternatively, closer to the axis, when y ^ 1 and xy <^1, 
we obtain the Fade approximant 

, , 1 — ,^ ^ , 

r{x,y)^— . (11) 



1 + xy 



This approaches the flat space limit as xy 0: 



Kij ~ Exp 



277r^ 1 
8 nfqf 



(12) 



Another regime of interest is at low energy, when y 
into the anthropic shell. In this limit [12] 



1: that of the final transition, down 



r{x,y) 



:i+x' 



{n, - 1/2)- 
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Figure 2: A plot of the flat space tunneling rate ( |T2l) . as a function of angular coordinates in 
three dimensions. The tunneling rate is represented here as being proportional to the distance 
from the origin. 

such that the gravitational correction is vanishingly small. As long as this condition is met, we 
can then approximate the tunneling rate down into the anthropic shell by equation f|T2l) . 

The fiat space tunneling rate is shown graphically, as a function of angular coordinates, 
in figure [2l It can be seen that Kij is significantly larger near the axes of this space, even 
in these few dimensions. In more dimensions, the contrast between maximum and minimum 
is increased. This expression for Kij will be used below to show that even a single transition 
(the last one) can produce considerable clustering of probability around the axes of the flux 
space. As the radius of the sphere being considered (and hence the energy level) is increased, 
gravitational corrections must increasingly be taken into consideration. The effect of gravity 
is to reduce the contrast between maxima and minima, smoothing the distribution of Kij and 
producing a plot that looks more like a sphere. 

As well as nucleating branes of single charge, one may also wish to consider the possibility 
of nucleating stacks of branes simultaneously. A study of such events has been carried out by 
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Garriga and Megevand [15]. In this work it is concluded that the instantons for nucleating 
stacks of branes may well not exist, but if they do then the tunneling rates for such events will 
be approximately given by B^j^ ~ kBij, where k is the number of branes in the stack and B^^'' is 
the exponent of the tunneling rate for nucleating the stack. This formula was found in the limit 
that the bubble wall is much less than the Hubble rate of the parent vacuum, in which case 
the thin-wall approximation can be used. When a sufficiently large stack is nucleated (again, if 
this is possible), then the thin wall approximation is violated and the bound on the minimum 
possible tunneling rate [L6\ is saturated. In this limit, the exponent of the tunneling rate was 
found to be B^^^^ ~ This rate is independent of the everything except the energy of the 
parent vacuum, and so once it is met all unlikely tunneling events proceed with this minimum 
probability, if instantons for these transitions exist at all. 

The possibility of nucleating stacks of branes significantly alters the dynamics of infiation, 
if it is allowed to occur. We consider it beyond the scope of the present study to perform a 
detailed analysis of this phenomenon here. Instead, we will first proceed using the standard lore 
of allowing only transitions between nearest neighbors in the BP landscape. After analyzing 
this situation we will then refiect on how nucleating multiple branes simultaneously will change 
this picture. We leave it to future studies to determine how likely these events really are. 

2.3 Anthropic Considerations 

In volume weighted measures of infiation the majority of the 3-volume of the universe, and the 
majority of new bubbles formed, are at very high energies. However, no observers can exist in 
these vacua. Firstly, the famous anthropic bound on the cosmological constant tells us that 
anthropically acceptable vacua must lie within the range < A < 10^^^° in order for structure 
to be able to form [17] (we neglect the possibility of small negative cosmological constant [18], 
for simplicity). Secondly, these high energy vacua decay far too quickly for structure to form. 
If we are interested in the likely locations in the landscape for observers to find themselves we 
should therefore not ask simply what are the most likely vacua, but instead what are the most 
likely vacua that (a) have a suitably low cosmological constant, and (b) are suitably long-liveci^. 

2.3.1 Cosmological constant bound 

As we mentioned in the preceding subsection, the condition (a) is tantamount to considering 
only vacua within a thin shell of width 10~^^°, and radius ^/2Aq. Any points further from the 
origin than this have A > 10~^^°, and so do not allow for the formation of structure. Any points 

^There are a number of other anthropic bounds that should also be applied from, for example, the fine 
structure constant and Newton's constant. For simplicity, and due to the limited nature of our model, we will 
not consider all of these bounds here. 
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closer to the origin are collapsing AdS vacua. The volume, V, of this shell is then given by 



V 



-l)\J 



The density of points in the flux space is (H/^'i) ^5 so that the condition for more than one 
point in this anthropic shell is 

2.3.2 Life-time bound 

The condition (b) requires that the time scale Tij = is less than the time required for 
structure to form, and life to evolve. This is clearly a very difficult time to specify precisely. 
For definiteness we will require the bound that the vacuum should have a lifetime > lOGyr, 
corresponding to 

The location of this boundary could be shifted, depending on the time at which observers are 
deemed to appear in any particular theory. Transitions out of any point that satisfies condition 
(a) will be adequately described by the fiat space rate, ( JT2l) . This rate is at a minimum along 
the 'diagonal', where Jnj^qf = 2Aq. The condition for suitable stability of these points is given 
by 

g.J^/^ > 13Af . (14) 

This condition on the slowest transitions rates is the requirement that there exist any points 
that live longer than lOGyr. The condition that all points in the anthropic shell are this stable 
requires a condition on the fastest transition rate, at the poles. Here nj^qf = 2Aq, and the 
tunneling rate ( fT2l) gives 

q. > 13Af . (15) 

The anthropic conditions (a) and (b) therefore give us acceptable ranges of the parameters 
qi and Aq, in terms of the number of dimensions J. Strictly, these relations should not be 
considered independent of each other as we require not only a minimum of one point in the 
anthropic shell and one region in this shell where points are stable, but one stable point in 
the shell. The combination of conditions (1131) and (1141) can then be considered a necessary 
condition, and the combination of conditions (IT^ and (IT3]) a sufficient one. These conditions 
are summarized in the table below, where it is assumed that all charges are approximately equal 
(so that all ^ g). 

It can be immediately seen, for reasonable parameters, that the condition for all points in 
the anthropic shell to be stable is often not satisfied. This means that the regions near the 
poles of this shell are not long-lived enough for observers to evolve. It can also be seen that for 
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some points 


all points 


enough points 




stable 


stable 


in shell 


J 


.Ao^^^ > 


.Ao > 


rX'' > 


100 


0.013 


13 


28 


250 


0.0033 


13 


8.3 


500 


0.0012 


13 


6.7 


1000 


0.00041 


13 


7.2 


2000 


0.00015 


13 


8.8 



fixed Ao (n!3|) provides an upper bound on q, and f|T4l) provides a lower bound. Similarly, for 
fixed q (fT3l) provides a lower bound on Aq, whilst (IT^ provides an upper bound. 

We expect the bare cosmological constant to be of order unity, and the charges g to be a 
couple of orders of magnitude below this scale. Reasonable values for the dimensionality of the 
phase space are 100 to 1000. For definiteness in what follows, we will often consider a model 
in which q = 10~^, Aq = 1/2 and J = 500. This satisfies the conditions ( |T3l) and (HM . but not 
( IT5|) . The condition (fTS!) could be satisfied by considering a smaller Aq. 



3 Populating at High Energies 

One of the benefits of volume weighted measures of infiation is that there exists a stable late- 
time distribution, which can be reached from a broad range of initial conditions: Once this 
'stationary' regime is achieved, the initial conditions of the system are quickly forgotten. All 
that is required for the universe to achieve stationarity is that at least one part of it should 
reach the highest energy states. Although this may require unlikely upwards tunneling events, 
the rapid expansion of the new higher-energy vacua produced can often be more than enough 
to offset the improbability of such events occurring. 

Of course there exists an upper boundary beyond which the semi-classical methods we em- 
ploy cannot be usefully applied. This is the Planck boundary, beyond which a full quantum 
theory of gravity is required for an adequate description of the space-time (if this concept 
survives). We will not consider such regimes here, but will instead restrict ourselves to consid- 
ering only sub-Planckian regions. We will now argue that on high energy shells, of constant A, 
probability will be distributed approximately isotropically. 

In a landscape with many high energy vacua the evolution of the universe, according to 
volume weighted measures, will not be dominated by a single high-energy vacuum, as may 
initially be suspected, but rather by the combined expansion of the many highest energy vacuslj. 

^We thank Andrei Linde for making this point clear to us. 



11 



The effect of the expansion of these vacua will give an exponentially growing probability 



which will be shared by the lower energy vacua in the manner outlined in section 2. Here Hh 
is the effective Hubble rate of the highest energy vacua, which we expect to be near Planckian. 
We may now ask ourselves what the condition is for a particular vacuum to find itself in some 
'smeared out' region at high energy, where all P ^ P/^. To answer this consider the evolution 
equation 

Pi — ^ ^ ^jiPj ^ ^ f^ikPi ~\~ "iHiPi 
j k 

~ ^hP. ^ t^jiPh - J2 ^'^^^^ + 3^*^- 

j k 

In going from the first to second line we have assumed that all neighboring vacua to Pi have 
approximately Pj ~ Ph-, i.e. that our chosen vacuum is in the smeared out region. This 
assumption will be consistent only if Pj ~ Ph, in which case the vacuum can consistently be 
described as being in a 'smeared out' region at high-energy. To see if this is true we can 
rearrange the equation above to obtain 

Ph 3AH + Z^^^l 

* 

where AH = Hh - Hi. Now, if 3 AH < ^ k^i, then 

In this case we have Pj ^ Ph, and the vacuum under consideration is in the smeared out region. 
(The order of magnitude here is due to the small difference in energy levels between neighboring 
high-energy vacua in the BP landscape). Alternatively, when 3 AH ^ X^^^i) have 

Ph 

— — ^ 1 

Pi E f^n 

and Pj ^ Ph- In this case the vacuum we are considering cannot consistently be assumed to 
be in the smeared out region. 

Furthermore, transitions between the highest energy vacua in the BP landscape will be only 
very weakly suppressed, with ~ 0(1). As such, the value of Kij/3Hh will also be 0(1) and 
the ratio of probabilities given by ([2]) will be 

^-0(1). 
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The products of any such decay will then be (approximately) as equally abundant as the vacua 
from which they decayed. If the highest energy vacua decay at these rates then the whole of the 
high-energy part of the landscape will very quickly become populated with an approximately 
flat distribution of probability. 

We have identified here reason for thinking that the highest energy levels of the BP landscape 
should have a probabihty distribution that is 'smeared out' (i.e. isotropically distributed). This 
smearing process erases any anisotropy, or favoring of any particular point or energy level in the 
initial conditions of the universe - as long as the highest energy vacua can be reached, we are 
left with a landscape that is populated at high energies by an approximately flat distribution 
of probability. Clearly this is a very different regime to the low energy one, in which tunneling 
events are suppressed by large factors of SHh/Kij ^ 1: Hence, there exists a boundary above 
which we can expect tunneling rapid enough to produce a flat probability distribution, and 
below which tunneling is heavily suppressed. 

4 Tunneling from High to Low Energies 

We will now consider how tunneling from high energies down to lower energies will proceed. 
First it is necessary to determine which paths will contribute the highest probability flux onto 
the anthropic shell. The shortest path from any point on the anthropic shell to the 'smeared 
out' surface, described above, will be a straight line extending in one dimension of the flux 
space. Any deviations from this shortest path will add extra factors of Kij/3Hh, which will 
suppress their relative contribution to the overall flux arriving at our point. It is now necessary 
to determine whether this suppression will be sufficient to beat the combined contribution of 
the potentially large number of paths which are slightly longer, and less favorable. 

The degree to which a longer path can be considered less favorable depends upon where it 
deviates from the straight line, shortest path, and whereabouts in the flux space the shortest 
path is located. For example, a straight-line shortest path which is parallel, and close to an axis 
of the flux space will have very little contribution from any other (longer) paths. This is due 
to the vacua along such a trajectory having one large flux, and the rest very small. Transitions 
in directions with small fluxes associated with them will then only progress at very low rates, 
compared to those in the direction of the single large flux. The shortest path, which is oriented 
in the direction of the large flux, will then consist of a series of rapid transitions. Any longer 
path will necessarily include additional transitions in a direction associated with a smaller flux, 
and so will be heavily suppressed. Away from the axes extra transitions at or near the surface 
beyond which rates are no longer suppressed will not incur a large penalty. We will treat this 
high energy surface as an initial boundary from which tunneling down to the anthropic shell 
begins. Any unsuppressed transitions before this surface can then be considered as contributing 
to the boundary, but not favoring any particular path thereafter. We will show that our results 
are not heavily dependant on the exact location of such a boundary. Extra deviations away 
from the shortest path at low energies will incur a large penalty, which will be enough to make 
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Figure 3: A schematic of the shortest path from high energies to the penuhimate sphere, parallel 
to the z axis. 



any such contribution to the total probability negligible. We therefore expect the shortest path 
from the high-energy boundary to be a good approximation to the probability of arriving at an 
acceptable point in the anthropic shell, both near and far from the axes of the phase space. 

We now need the appropriate transition rates along these shortest paths. Although the 
approximation to the gravitational correction ([9]) is the appropriate one for most of the high- 
energy transitions (in a high dimensional landscape), it is not the one in which we are interested 
here. A straight-line extending on one dimension of the flux space, from the anthropic shell to 
high energies, will always be contained within a cylinder about one of the axes, of cross-sectional 
area less than that of the anthropic shell itself. The appropriate gravitational correction is 
therefore given by flTT]) along such a path, until the last transition where it is given by f[T^ . In 
this section we will consider the regime where (fTT!) is valid. The effect of the last transition will 
be considered separately in the next section. 

The tunneling rate prescribed by ffTTj) is 



and the total probability of tunneling down from high energies to the penultimate sphere, before 
the anthropic one, is then given by the product of these tunneling rates along the shortest path, 
shown schematically in figure [3l 

The value of Aj along this path is 




A 1 2 2 -'-2 
Ai = iz^i^i - 7;Zo, 



(16) 
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where Zq = 2Aq cos^ and is shown in figure [3l The value of riiqi, for any i, is 



UiQi = qi{noi +j) = zo + jqu 



(17) 



where noi is the value of Ui in the anthropic shell, and j is the number of fiux quanta in z away 
from that point. Now the total supression along this path is given by 



n 



Hi 



Exp 



N 



4n,^g/ - 3A, 
Anfqf + 3Ai J n^q^ 



where Aj and nj^j are given by fll6p and fll7p . and N is the number of transitions involved in 
tunneling from high energy to low energy. In a continuum approximation this expression can 
be approximated by 



kt — Exp 



Exp 



27n\ Z-^+i/^ 
8 J3/2 



4:nfqf-3Ai \ dj 
Anfqf + 3Ai 
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3(nfg2_2A,) 2njq^ 
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As we have already discussed, the leading contribution to this quantity will be the low energy 
limit, where extra transitions are heavily suppressed. The contribution from extra transitions 
at high energy will be negligible, and correspondingly we take ^ 00. In this limit we then 
have 



— Exp 



3^2 



In 1 



I2zl 



.0 V ll(2zo + 3gi)V (2;2o + 3gi)VJ 
The limit A^ ^ cxo is further justified in appendix |Al where integrations are performed to finite 
values of A^. It is found that the quantities we are interested in are largely insensitive to the 
exact location of the high-energy boundary, and so we are justified in approximating it to be 
at infinity. For zq ^ qi the expression above becomes 



~ Exp 



277r2 /8 



Exp 
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'0 J 



(19) 



It can be clearly seen that the tunneling rate is significantly faster to points in the anthropic 
shell with large z (i.e. close to the axes, away from the diagonal). This is due to the supression 
associated with the extra transitions required to get to the points further from the axes. We 
will now use this expression for the tunneling rates to determine where on the sphere a 'typical' 
vacuum is located. In order to do this we will need to balance the faster transition rates at the 
poles of the sphere, with the increased volume along the diagonals. 
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4.1 Expectation Values 



In order to consider the distribution of probability on a sphere it is useful to determine the 
moments 

(20) 



if) ^ J fnrdn 



where / = f{xi, z) is any function of Xi and z on the sphere, and dVL is the surface element of 
a hypersphere. This expression takes into account not only the effect of the transition rate kt, 
but also the effect of increasing surface area away from the axes. 

We can consider the penultimate sphere as being divided into 2 J symmetric sectors, each 
centered around one of the 2 J axes of the space. Each of these sectors corresponds to the 
region a on the surface of the penultimate sphere which is closest to the axis about which it is 
centered, z is then given, in each respective region, as the corresponding value of the Cartesian 
coordinate along its central axis, and the whole sphere is built up by evaluating {-n?) in any one 
of the regions and multiplying it by 2 J. The integral that we are now required to evaluate is 



1 j-i 2 /j-i \ 

J dzHj dx,fe~''/^'s\y2^' + ^'-n^ 



(21) 



where we have taken Aq = 1/2. The Cartesian coordinates Xi correspond to all of those 
directions which are orthonormal to the central axis z. The limits of integration, —z and z, 
then pick out the sector discussed above. Included in the integrand is a delta function intended 
to select the surface at ~ V2Ao = 1, which corresponds to the anthropic shell. Here (p) is 
being evaluated in any one of the 2J symmetric regions. We will be interested in the choices 
f = z and f = X, where x is any Cartesian coordinate orthonormal to z. The moments (z^) and 
{x^) then give us the expectation values of these quantities, and their ratio gives us a measure 
of the peaking of probability on the anthropic sphere due to tunneling down from high energy 
to the penultimate sphere before the anthropic one. 

In appendix [B] we find that the expectation values for (z'^) and (x^) can be given in terms 
of one-dimensional integrals as 



r(^) .0 



dyy"^ 



^E(v/i2T7) 



J-1 



,-(12+!/2) 



:221 



and 



n^j JO 



dyy"^ 



2(12 + 2/2 



-(12W) 



X 



12 + y^E{^VlT^) 



12 + y2 



J-1 



(23) 
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Here the function E{p) is the integral of a Gaussian distribution, or error function, given by 

2n-y^j;e-'^'dq. 

The measure of peaking that we are interested in is given by the ratio of these two ex- 
pressions, {z^)/{x^). This is shown graphically in figure HI It can be seen from this plot that 
the peaking produced by tunneling down to the penultimate sphere is extremely sharp, and 
increases with the dimensionality of the flux space. A fiat distribution of probability would 
produce the result (z^) / (x^) = 1. The large values of this quantity here show that there exists 
a strong peak in probability on the anthropic shell. 



4.2 The Single Peak Approximation 

We are now in a good position to investigate approximations that will make these quantities 
more transparent. A particularly useful simplification is to extend the domain of integration of 
equation fl2Tl) . to make it independent of z. If we were to allow the limits of integration for the 
coordinates Xi to run from —1 to 1, instead of — z to z, then this would allow us to transform 
to a hyper-cylindrical coordinate system in which J — 1 of the dimensions could be integrated 
over in a straightforward fashion. This would leave a simple one-dimensional integral that can 
be solved by saddle point approximation. We will call this the single peak approximation, as 
we will be integrating over a larger region that should contain several peaks of probability; we 
will, however, ignore all but the single central peak of the initial region and show that this is a 
good approximation to integrating over the correct domain, as in fl21l) . The reason why this is 
a plausible approximation to consider is that the rate of transitions f|T9|) can be seen to become 
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very small at low z. The extra region being integrated over in the single peak approximation 
is all at low and so it is a reasonable proposition to ask if integrating over such an extended 
domain yields similar results to integrating over the actual domain of interest. We will argue 
that it does. 

The single peak approximation yields the results 

"OO J— 1 /.oo 







dz^ 



dzz 



oo 



Vl2 + ^2 



J-1 



(24) 



which can be compared to equation ( !22l) . and 

J-1 



' J+2^ 



' J+2^ 
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dxiX^e'^^^ x2+22)(i+i22-2) 



2(12 + ^2) 



V12 + z- 



:E(V12+^) 
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,-(12+^2) 



(25) 



which can be compared to equation fl23l) . We denote quantities evaluated in the single peak 
approximation by an over-bar. The fractional error of considering this approximation can then 
be written as 

_ {z')/{x') - {z')/{x') 

= — wim — ■ ^ ^ 

This quantity is vanishingly small for all J, showing that the single peak approximation is very 
good. 



4.3 Analytic Approximations 

Using the single peak approximation we will now find analytic expressions for these moments. 
Recall the tunneling rate is given by 



kt ~ Exp 



12" 


= Exp 


K' 










7^. 



(27) 



where K = 12, and to be concise we have replaced zq with z. The normalization of the 
probability distribution is then 



/„i J-1 „i / J-1 

Kdn = 2 J y n y dxie-^^''5 ij^x'^i+z^-l 



{2i 
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In the single peak approximation, the hmits on the Xi integrals are replaced with -1 and 1. 
Changing coordinates from Cartesian to hyper-cylindrical, and defining = 
rewrite as 



N = 2JVj^2 [ dz [ drr^-^e-^/''6{r^ + - 1) , (29) 
Jo Jo 

where Vj_2 is the volume of a ( J — 2)-sphere. Eliminating the delta function with the r-integral 
leaves us with an integral over z: 

N = 2JVj.2 f dzil - z^)'^e-^'^' ~ 2JVj^2 f dze<'"''"~i'^^^-'"^) . (30) 
Jo Jo 

This integral can be approximated by a saddle point, which is at 

At first glance it may seem that this result is in contradiction with the validity of the single 
peak approximation. After all, z is small for any reasonable (C < 1) and therefore r is close 
to 1. However, although z is small, a typical Xi (~ l/\/~J) is much smaller and so the peak in 
the distribution on the sphere is, in general, closer to the 2;-axis than to any of the other Xi 
axes. This issue is discussed in more detail in section |6l 

We can similarly calculate the expectation values of powers of the fluxes: 

N Jq N Jq 

This should be normalized using the A^ calculated above. The integrand is again extremized at 
z ~ (^)^^^ (with an 0{{2K / jy/^{k/ J)) correction). Substituting this back into the integral 
above, and evaluating the saddle point, gives us 

,,,j^4r'-inT ■ (33) 



J J \J J 

Here the normalization A^ cancels the exponential, the remaining integral and various other 
numerical factors. In a flat measure the same expectation value is given by 



{z%.t = ^ / dr{l - r'Y'^-^yV-^ . (34) 

^Vflat Jo 



oo 



The scaling with J can again be obtained, and is now found to be 



~ ( J ) . (35) 
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Symmetry considerations then give the expected values of the remaining fluxes as 



i^') - -jm^ (36) 



and a typical flux vector can be estimated as 




J 'jT7^'---'ji7^)' (37) 

which, for J = 500 and q — 10~^, is 

n ^ (45 , 4.5 , . . . , 4.5) . (38) 

Therefore, the most likely vacua produced by tunneling down from high energies to the penulti- 
mate sphere arc those with one flux approximately ten times as large as the others. The spread 
of probability about this most likely point will be addressed in section 6. Here we have assumed 
that all charges are the same. Allowing different charges will produce sharper peaking near the 
axes associated with smaller q. 



5 Transitions into the Anthropic Shell 

Having calculated the probability of transitioning down onto the penultimate sphere, above the 
anthropic one, we will now evaluate the effect of the last jump. The flnal transition rates into 
the anthropic shell are well approximated by the flat space rates 

KF ~ e-^''' (39) 

where z denotes the value of the Cartesian coordinate in the direction of transition, and C is a 
constant {C ~ 0.3, when q^ = 10^^). Strictly, z should be a discrete number and correspond to 
the penultimate shell, one step up from the anthropic one. In our continuum approximation, 
however, z is a continuous variable, which we can evaluate on the surface of the anthropic shell. 

Transitions will occur in all J dimensions of the space. However, we see that Kp is very 
strongly dependent on the value of z, allowing us to approximate the probability density on the 
anthropic sphere as being primarily due to transitions from the largest of these z. We again 
consider the sphere as being divided into 2J symmetric sectors. The integral that we are now 
required to evaluate is 

if) = I' dz n £ dx,fe-^/^'6 + - ij- (40) 

The Cartesian coordinates Xi correspond to all of those which are orthonormal to the central 
axis z. The limits of integration, —z and z pick out the sector of interest. 
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In appendix Owe show that the single peak approximation for the final transition is accurate 
to better than 0.03% at J = 500, and that this accuracy increases as J is increased. Having 
justified the validity of this approximation, we will now use it to obtain analytic expressions 
for the peaking of probability produced by the final transition. 



5.1 Peaking From the Last Step 

Working in the single peak approximation we can use saddle points methods, as before, to 
estimate the level of peaking from the final transition. We begin with the integral 

/„1 J-l „2 / ,7-1 \ 

Kdn = 2Jj dzYlj dx^e~^^''5 l^x^i + -ij , (41) 

which has its integrand extremized at 

Again, the integral has most of it's support relatively close to the z axis, consistent with the 
single peak approximation. 

Once again we can calculate the expected fiuxes from the integral 

(^fc) = %i dr{l - ^2)(fc-i)/2^J-2^-c/(i-.^)3/^ ^ (43) 







which has its saddle point given by|42] (with an 0{ek'^/J) correction), and thus 

'^/S . fc/5 



(44) 



This allows us to write a typical vector as 



~ (jTTF ' JTT^ ' • • • ' 7171) ■ (45) 

For J = 500 and q = 10^^ we should then typically expect to find one set of preferred fiuxes to 
be over 6 times greater than the others: 

n~ (29 , 4.5 , ... , 4.5) . (46) 

These saddle point approximations are shown graphically in figure [5], where they are compared 
to direct numerical integrations. 
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Figure 5: A plot of calculated numerically (solid line), and by saddle point approxi- 

mation (dashed line). The saddle point approximation is taken from fHSl) . and can be seen to 
be a reasonable approximation to the numerics at large J. The accuracy of the saddle point 
approximation increases with J. 



The peaking found here is clearly less than that found in the previous section, and we will 
now show that its effects are negligible in comparison. The total tunneling rate from combining 
([T9D and ([39]) is 

12 0.3" 



K Exp 



^0 ^0 



It can be seen that the second term in this expression only becomes important at about zq ~ 
1/40. Before this the first term dominates, and the expression is well approximated by (fT9|) . 
This critical value of zq can be seen to be an order of magnitude smaller than the location 
of the saddle points involved in either this, or the previous, section. As both of the saddle- 
points are located comfortably in the region where ( JT9l) is a good approximation to the total 
tunneling rate, we consider the effects of this last transition to be negligible in comparison to 
those considered in the last section. This is supported by the relatively mild peaking found 
above. 



6 The Most Probable Vacua 

Having demonstrated that regions of the anthropic shell are prefentially selected by volume- 
weighted measures, we now turn our attention to how sharp this peaking is. We begin by finding 
analytic estimations for the width of the probability distributions, and then use these to obtain 
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characteristic numbers for the volume of the universe that exists in some small fraction of the 
anthropic sphere. We then consider the issue of life-times. By imposing an anthropic limit on 
the minimum life-time of acceptable vacua we find that there exists a pressure favoring shorter 
life-times. 



6.1 Peaking About the Most Likely Point 

We will now refine our assertion that the measure on the anthropic sphere is sharply peaked, 
and give a demonstration of the level of this peaking. We begin our argument by considering 
the tunneling down from high to low energies, as discussed in section HI We do this by noting 
that near the peak the distribution is Gaussian, and then calculating the standard deviation. 

We begin with the integral (I30p and note that the saddle point approximation gives us not 
only the position of the peak, but also the width. The leading contribution to this is found to 
be 

Now consider the fiat probability distribution on the sphere, where 

iVfiat = K7-2 / dz{l - ^2)(^-3)/2^ (48) 



at 2; <^ 1 this can be approximated by 

iVflat ~ Vj.2 j (i^Exp[-JzV2] . (49) 

Evidentally the integrand falls at least exponentially for z > J^^, and grows rapidly above 
z < . Roughly speaking, this means that whilst the probability distribution integral has 
most of its support around z'^ = 24/ J, most of the volume of the sphere lies further out, around 
z = 1/J. This is represented graphically in figure O where it is clear that the overlap between 
the probability peak and the volume peak decreases as J is increased - i.e. more and more 
of the volume of the anthropically acceptable parts of the multiverse come from smaller and 
smaller fractions of the surface of the anthropic sphere. 

We can then estimate what fraction, /, of the sphere is covered by mcr of the probability 
distribution, by evaluating the volume integral (jUj) between the limits {z) — ma and (z) +ma: 



(2>+m(8J)-i/2 



/ = 2J-^ / d^Exp[-Jz72] . (50) 

We can bound this integral by multiplying the maximum of volume integrand (which is at z = 
(2Ky^^ / J^^^ + m(8J)~^''^) by the total range of the integral (2ma). The leading contribution 
is given by 

J ~ 2J— — e ' ~ 2Jy7 — e 2 ~ Amy Je 2 . (51) 

Vj_i J Vj-i J 
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Figure 6: A sketch of the peaking of the probabihty measure as a function of 2;, along with the 
volume of the sphere. The overlap between the probability peak and the volume peak decreases 
as J is increased. 



For J = 500 we find that 99.9999% of the volume of the universe (5 standard deviations) covers 
~ 10^^^ of the anthropic region. 

This argument can be repeated for the final transition down onto the anthropic sphere, 
discussed in section |5l Here the fraction of vacua on the anthropic sphere within ma of the 
peak is 

Vj 9 2m .73/5 ./^ .73/5 

/ ~ 2J^^j^e-^ ~ 4mj'/'e-^, (52) 

where a ~ we then find that 99.9999% of the volume covers only ~ 10 ^ 

of the anthropic sphere. Evidentally, this is considerably less than that due to tunneling down 
from higher energies. 



6.2 Pressure Against the Life-Time Boundary 

In the preceding sections we have found the most probable vacua that have a cosmological 
constant in the anthropically acceptable range < A < 10^^^'^. However, for observers to 
appear we require not only an anthropically acceptable value of A, but also a vacuum that is 
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long enough lived for life to evolve. In section 2.3 this was taken to be ~ 10 Gyr. The vacua 
located at the peaks of probability (found in the previous sections) do not satisfy the bound 
r > 10 Gyr, and so are unlikely to be stable enough for observers of our type to occur. If 
we wish to determine the probability of vacua that have both an acceptably small A, and an 
acceptably large life-time, then we must calculate probabilities only up to the boundary at 
which T ~ lOGyr. Using our knowledge that there exists a sharp peak of probability outside of 
this region, we should expect to find a selection pressure pushing against the anthropic life-time 
barrier (so that vacua closer to the peak are favored). Here we will investigate how strong this 
selection pressure is, and use it to make a prediction for the half-hfe of acceptable vacua. 

In order to establish the extent to which vacua near the life-time boundary are favored we 
will consider a thin region next to it. The life-time boundary is at a constant value of z, and 
we will label it by Zi. The thin region we are considering next to this boundary will extend 
down to a constant z = Z2, where Z2 < Zi. The volume contained within this region can be 
approximated by 



where Vj_2 is the volume of a (J — 2)-sphere. The probability of finding a vacuum in this region 
can then be estimated using the formula (fT9|) . and by integrating over the thin band: 



These integrals can be performed numerically, the results of which are shown in figure [71 Here 
we have normalized AP so that AP/P ^ 1 as 2:2 ^ 0. The volume AV is presented as a 
fraction of the volume of the unit sphere. It can be seen from the figure that the probability is 
pressed up against the boundary, within a fractional volume of ~ 0.001. 

In order to determine the typical distance of an anthropically acceptable vacuum from this 
boundary we will find the width of strip that is required for half of the probability to be within 
it, and half outside of it: The half-life of the distribution. The value of Z2 for such a band is 
~ 0.99952;i. The life-time of vacua at this value of z is ~ 12.2 Gyr. Given the model we have 
been investigating, with the parameters chosen and assumptions made, we should therefore 
expect a typical anthropically acceptable vacuum to decay after about 2.2 billion years from 
the moment when the first observers appear (at lOGyr). 

One may have initially suspected that the exponential pressure up against the boundary 
should cause decay of vacua after a very short time. However, the exponential sensitivity of 
life-time on z counters this pressure, and results in vacua being reasonably long lived. 

7 Discussion 

In this paper we have considered the application of volume weighted measures of eternal inflation 
to the Bousso-Polchinski landscape. We have found that the cosmological dynamics of volume- 
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Figure 7: A plot of probability, AP, in a thin band of volume ISV next to the life-time boundary. 
Here we considered the model J = 500, q = 10~^ and Aq = 1/2. 

weighting results in significant peaking of probability on the surface of the anthropic shell. The 
anisotropy of tunneling rates in these models leads to the selection of vacua that have one flux 
approximately ten times larger than the other 499 (for J = 500, q = 10~^ and Aq = 1/2). 
The peaking of probability about this most likely outcome is very sharp, with 99.9999% of the 
probability on the anthropic shell expected to be in a region that is only ~ 10^^^ of its surface. 
The preferential selection of a single large flux might suggest that warped compactifications are 
favored by this measure. This is in contrast to the results of Bousso and Yang [9] who find that 
the holographic measure [8] selects fluxes along the diagonal of the flux space, with no one flux 
being especially large. 

Although this peaking gives us a sharp prediction for the flux configuration we should expect 
from volume-weighted inflation, it still allows for a satisfactory explanation of the cosmological 
constant problem. The reason for this is the very large number of points that are typically in 
the anthropic shell, ~ iQ^®^ for the parameters chosen above. Even an area as small as 10~^^ of 
the anthropic shell still then contains ~ 10^^^ vacua with anthropically acceptable values of the 
cosmological constant^. This is more than enough to produce the desired smooth distribution 
in A that is required for an acceptable anthropic explanation of the problem. We do not expect 
the volume weighted measures of inflation to produce a staggered distribution, such as that 
found in the bubble counting measure of [B] for some values of parameters [7], [T^ . 

We have also found that the most likely vacua in the anthropic shell often have a life-time 
too short for observers, such as ourselves, to appear. This additional anthropic pressure makes 

^Other reasonable values of parameters also produce acceptably large numbers. 
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the vacua that are sufficiently long-hved, and as close as possible to the favored points, the 
most likely for observers to find themselves in. Such a configuration selects vacua that are 
long-lived enough for observers to occur, but short-lived enough to be as close as possible to 
the overall peak of probability. We find that if it takes 10 billion years for observers to occur, 
then the most likely vacua for them to find themselves in will decay a further 2.2 billion years 
after this time, in our simple model. This result is dependent on the decay rates in section 
2.2 being applicable to every vacuum in the anthropic shell. If anomalous vacua with longer 
life-times (due to e.g. super-symmetry) are identified, then these vacua would be immune to 
this selection pressure and may well prove to be more likely, depending on their number and 
flux configuration. 

This apparent prediction of a doomsday might provide a plausible explanation of the "co- 
incidence problem" of cosmology. This problem, loosely stated, asks why it should be that 
the energy densities observed in dark energy and matter should be comparable at the present 
time. If our pocket of the universe is allowed to exist forever, then a randomly distributed 
observer should expect to find themselves at a time when the dark energy dominates, and all 
other energy densities are negligible. The existence of a doomsday, a few billion years after 
the ffist observers appear, would change this picture and make it considerably more likely for 
a random observer to find comparable values of the dark energy and matter densities in their 
observable universe. A number of other solutions to this problem have also been proposed in 
the literature. For example, a compelling and model independent solution exists in the context 
of holography |20j . 

The results above have been calculated under the assumption that only tunneling events 
to nearest neighbors in the landscape are possible. However, it may be the case that multiple 
branes can be nucleated simultaneously, in which case we would have to be re-evaluate our 
results. For example, according to the results of Garriga and Megevand [TS] the probability 
for simultaneously nucleating a stack of m branes, if possible, is approximately given by k"^, 
where k is the rate of nucleation of single branes. If a result of this kind can reliably be applied 
to the situation we have been considering here, then it will have important consequences for 
our results - such an event, although unlikely, would be more likely than nucleating a single 
brane, and then m — 1 subsequent further branes from the lower energy vacua produced. To 
see this notice that the tunneling rate flTIJl) . for transitioning from the highest energy vacua 
to low energies, would become kt ~ g-is/mzg gg^-g -yyg have multiplied the numerator of the 
exponent by m, as suggested by [15], and multiplied the denominator by m^, as the larger 
jumps produced by nucleating multiple branes are equivalent to rescaling zq uizq. A large 
m can now be seen to decrease the exponent of kt-, therefore reducing the peaking. Clearly 
the possibility of multiple brane nucleation events, if allowed to occur, will have a potentially 
important effect. 

Multiple brane nucleation events could also have a second important effect, related to the 
minimum tunneling rate alluded to in section 2.2. If it should be the case that any number of 
branes can be nucleated simultaneously, then even the most unlikely events cannot occur at a 
rate less than e~^'^ if they are allowed to occur at all. For vacua with sufficiently large A 



27 



this rate is not negligible, and tunneling events to distant corners of the landscape could occur 
with a frequency comparable to that associated with vacua that are much closer. Clearly, in 
such a regime the sequences of tunneling events discussed above make very little sense; instead 
one would expect a very nearly fiat distribution of probability over the entire landscape as the 
highest energy vacua decay to all other vacua at approximately the same minimum rate. This 
is a very different picture to the one we have been focusing up in this paper. If only one brane 
is allowed to be nucleated at a time we have very strong peaking of probability on the anthropic 
sphere, if, however, any number of branes can be nucleated together then we have no peaking 
whatsoever. As it is this peaking that is responsible for a preference for short-lived vacua, 
removing it also removes this pressure. In these situations it would be the combined effect of 
the many unlikely transitions that would be of greatest importance. 

Another simplification that we have employed in this study is to neglect the necessary 
period of slow-roll at the end of inflation, which is not included explicitly in the Bousso- 
Polchinski landscape. Linde argues in [5J that volume weighted probabilities should be weighted 
by the amount of volume produced during the slow roll period. Such exponentially large 
factors could swamp the effects we have discussed here. Clearly it will be important to have 
a proper understanding of how to include slow-roll in this picture, in order to gain a complete 
understanding of how volume-weighted inflation proceeds in a more realistic landscape. 
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A Changing the Location of the High-Energy Barrier 

In section 4 the tunneling rates to the penultimate sphere, above the anthropic one, were 
calculated in the limit of the high-energy 'smeared out' surface being at an infinite distance 
from the anthropic surface. From the results of section 3 it is clear that this surface should, in 
fact, be at a finite distance, and so we must justify our approximation of it being at infinity. 
To achieve this we will take some plausible values for the actual location of this boundary, as 
suggested by section 3, and use them to evaluate the quantities we are interested in. We will 
find that the results of our calculations are relatively insensitive to the precise location of the 
high-energy surface, so that we are justified in approximating it to be at infinity. 

Let us first consider a boundary defined by the condition 3 AH = Kii- This condition can 
be considered as an upper estimate for the high-energy surface: Vacua that meet this condition 
are still within the 'smeared out' region. Using the high energy tunneling rate flTT]) for n^i, this 
condition picks out a surface that can be used as an upper limit for the integral (fT8l) . This 
surface intersects the 2;-axis at z ~ 1.23, and drops to slightly lower z as you move away from 
the axis. 
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Alternatively, we may wish to consider a high-energy boundary defined by 3 AH = 10Kj|. 
Vacua that meet this condition have violated the condition to be in the smeared out region 
by an order of magnitude, and so can be considered to be comfortably outside of the smeared 
regime. This condition is then a lower estimate for the high-energy surface. Again, using ffTTl) 
this condition picks out a surface which can be used as an upper limit for the integral f|T8|) . 
This lower estimate cuts the anthropic shell, and intersects the z-axis at z ~ 0.61. Again, this 
surface drops to lower z as you move away from the axis. The intersection of this surface with 
the anthropic shell clearly gives a different picture - now tunneling to vacua in the anthropic 
shell, above this surface, can be tunneled to in an unsuppressed way. Vacua below this surface 
must be reached in the usual way, by a downward series of suppressed jumps. 

Using these two estimates for the high-energy surface we can numerically find the quantity 
we are most interested in: The saddle-point of the integral J n^idfl. This quantity determines 
the position of the peak of probability, and allows a straightforward estimate of the moments 
{z'^) and {x^) (see section 4.3). The location of the saddle-point in these estimates, together 
with the case of a high-energy boundary at infinity, are given in the table below. 



Location of 


low 


hi{ 


^h 


infinity 


high-energy surface 


{3AH = 10Kj,i) 


{3AH 




{z oo) 


Location of 










saddle-point 


zo = 0.433 


Zo = 


0.439 


Zo = 0.440 



It can be seen that the location of the saddle-point is only weakly dependent on the location 
of the high-energy surface. The difference between our higher estimate for the 'smeared out' 
surface, and taking this surface to be at infinity, is only at the level of ~ 0.2%. We attribute 
this to the fact that most of the suppression is due to transitions at low energy, so including 
extra transitions at high-energy gives only a modest contribution. The difference between our 
lower estimate, and taking the high-energy surface to infinity, is still relatively small at the 
level of ~ 1.6%. This is surprisingly accurate considering that in this approximation some of 
the transitions are completely unsuppressed. We attribute this accuracy to the unsuppressed 
points near the z-axis contributing only a small amount to the total probability fiux, due to 
the small amount of area near the pole compared to the large area away from it. Further from 
the pole, where the area is much larger, the high-energy surface is still far enough away to be 
well approximated by being at infinity. 
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B Expectation Values Due to Tunneling From High En- 
ergy 

To evaluate the expectation values {z^) and (x^), from equation (!2T|) . consider the integral 

/"OO J — 1 nZ 

1= dzT\ dxiiy2x^,+ z^)e-^^^''^+''^fKT (55) 

^0 , J-z 



oo 



,.2 







where / = f{xi,z) is any function of Xi and z on the surface of the sphere. In moving from the 
first to second line we have simply transformed to a set of hyper-spherical polar coordinates, 
where r = xf + z"^. In order to evaluate this second line, and obtain the third, we have had 
to be careful in how we have defined kt', specifically, we have required that kt = kt{^)- This 
may at first seem contrary to ( |T9l) . but it should be remembered that we are only interested in 
the value of on one particular surface - that of the anthropic shell. We can then re-write 
(fT9l) as 

where z = z{Q) is the value of z on the unit sphere at the point described by the polar 
coordinates Q. For points not contained on the unit sphere, kt is obtained by connecting the 
point to the origin with a radial vector and reading off the value of z as this vector crosses the 
surface of the unit sphere. Similar triangles then give 



or 



Similarly, we can re-write as 



1 P 

^^_e-i2(E,-?+^^)/^^ 

Y,i xj + z"^ x^ 



1 x^ 

where x = x{Q) is the corresponding value of x at the point a suitable radial vector crosses the 
anthropic shell. With these new expressions for z"^ and x^, the integral (I55p can be written as 

j-i 



c\ poo ^ pz 



Taking f = z the integrals over the coordinates Xi can all be performed exactly, giving the 
moment (z"^) as 



.2 



-I 



VI2T 



-.E{VT2 



Z' 



2) 



(56) 
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where E{p) is the integral of a Gaussian distribution, or error function, given by 27r~^/^ e~'^^dq. 
This analysis has reduced (z'^) from a J dimensional integral, which is potentially very difficult 
to solve, to a one dimensional integral that can be solved by saddle point approximation, or by 
a simple numerical calculation. 

We can proceed in a similar fashion to obtain an expression for (x^), which can also be 
expressed as a one-dimensional integral. Such an analysis gives 



J-1 



- 2 ^ 




^/7TZ 



_Vl2 + z^ 



E{Vl2 + z^) 



J-2 



,-(a;2+22)(l+l22-2) 



2(12 + z^) 7?V12 + z^E{Vl2 + z^ 

y/nz 



_Vl2 + z^ 



E{Vl2 + z^) 



j-i 



-(12 + 22) 



(57) 



Going from the first to the second line, we have integrated over all of the Xj dimensions, except 
for X (that is, the particular we are calculating the expectation value of). Going from the 
second to third lines we then integrate over x. 



C The Single Peak Approximation for the Final Transi- 
tion 

The integral fHOj) is not trivial to evaluate, and so we will briefiy consider a toy model in which 
the calculations are significantly simpler to perform. This model will allow us to fix ideas and 
explicitly investigate approximations which will be of direct use in solving ( l40l) . The toy model 
we propose is one in which the rate of transitions goes as 

ntoy e-^^/'\ (58) 

where (7 is a constant. This rate is clearly the same as (fT9ll . with a different numerical factor in 
the exponent. We already know such models can be solved exactly, from the preceding section. 
We will use this knowledge to justify the single peak approximation for the tunneling rate fl39p . 

In order to show that the single peak approximation is valid for fl5^ . we will proceed as 
follows. First we will investigate this approximation in our toy model. Having shown the 
extent to which this approximation is valid in our toy model, we will argue that it should be 
a significantly better approximation in a more realistic model. We will then show the results 
of numerical calculations of the more realistic model at low dimensions. This will support our 
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5 




Log J 



Figure 8: A plot of 5 as a function of dimensionality of the space, J, in the toy model (dashed 
line) and in a more realistic model (solid line). The value of C, in the toy model is taken to be 
0.7, and the value of C, in the more realistic model is 0.3. 



previous argument and, when extended to larger dimensions, should provide credible evidence 
that the single peak approximation is good. 

For our toy model, the single peak approximation yields the results as for 0191) . with 200 
replaced by C. We now plot S for the present case as a function of J in figure [HI It can be seen 
from the figure that the error, 6, can be as large as ~ 0.3 at low dimensionality {J 2). At 
larger dimensionality, however, the value of 6 can be seen to drop off, and for J = 500 can be 
seen to be less than 0.03%. We expect that this is due to the fact that at higher dimensionality 
the boundaries of the domain of integration of the actual region of interest extends to lower 
and lower values of z. This can be seen straight-forwardly from the expression for the diagonal 
(that is, the place furthest from all axes where the values of all coordinates are equal): 

J 1 
J]x^ = l ^ ""'^T/J' ^^^^ 

i ^ 

As this region extends to lower z, it covers more of the region that is integrated over in the single 
peak approximation, making the approximation increasingly more accurate as more dimensions 
are added. 

Having established that the single peak approximation in our toy model is accurate to better 
than 0.03%, when J = 500, we can begin to consider the accuracy of this approximation in 
more realistic models. We expect that the single peak approximation should be more accurate 
in the more realistic models. The reason for this is that the rate of transitions is a steeper 
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function of z {e~^l^^ instead of e^^/^^). With this in mind we should expect the extra area 
integrated over in the single peak approximation to make even less of a contribution, as this 
extra area is at low z. 

The calculation to be performed in the more realistic model is particularly cumbersome 
at large J, as we do not have the pleasant reduction to a one-dimensional integral that was 
exhibited in the toy model. For low-dimensionality, however, numerical evaluations can be 
performed. The recipe in this case is similar to that of the toy model. Consider the integral 

dzHj rfx,(^x2 + ^2)/V(^>^?+^')''\^ (60) 



oo 



Now, by the same method as before, we write the transition rate as a function of angular 
coordinates 



Ki? ~ e 



and 



so that fl60l) gives 



1 p 



Q poo ^ nZ 

(h = ^J dz]J dxJ^e-^^^^^-^'^'^'^^-<^^-^^ (61) 

i { — ) Jo i J-z 

In the single peak approximation this integral becomes 

Q roc J^l /*oo 

i {. — ) Jo j i-oo 

This last integral can be written in a significantly simpler way in terms of a set of hyper- 
cylindrical coordinates (which has been the reason for this whole subsection). Here, however, 
this simplification is not particularly beneficial as in order to evaluate 5 we must also analyze 
the more complicated integral (!6T!) . over the actual region of interest. 

The results of numerically integrating fl6Tl) and fl62l) are also shown in figure [HI Here 5 
is defined in the same way as in the toy model, and we have set C = 0.3. The value of C, 
for the toy model, was taken as 0.7 in order to match the two models at J = 2 (the lowest 
dimensionality for which 6 is well defined). It can be seen that the magnitude of 6 in this more 
realistic model is indeed smaller than in the toy model, as conjectured earlier. Although we 
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do not have reliable estimates of 5 at large J in the more realistic model (due to the difficulty 
of numerically integrating over a large number of dimensions), we argue by analogy that we 
should expect 6 to continue dropping as J is increased. 

This concludes our argument in favor of the validity of the single peak approximation. We 
have shown that there exists a toy model in which (/^) can be evaluated in a more straight- 
forward fashion. In this toy model we have shown that the single peak approximation is 
accurate to ~ 0.03% at J = 500, and that this accuracy increases as J is increased. We have 
then argued that in our more realistic model the single peak approximation should be more 
accurate than in the toy model. Numerical investigations at low dimensionality indicate that 
this expectation is correct, and comparison with the toy model leads us to expect that the 
single peak approximation in more realistic models should indeed be a good approximation at 
large J. 
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